Dynamics of an SO(5)-symmetric ladder model 
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We discuss properties of a recently proposed exactly 5'0(5)-symmetric ladder model. In the 
strong coupling limit we demonstrate how the SO(3)-symmetric description of spin ladders in terms 
of bond Bosons can be upgraded to an S'0(5)-symmetric bond-Boson model, which provides a 
particularly simple example for the concept of 50(5) symmetry. Based on this representation 
we show that antiferromagnetism on one hand and superconductivity on the other hand can be 
understood as condensation of either magnetic or charged Bosons into an RVB vacuum. We identify 
exact eigenstates of a finite cluster with general multiplets of the 5*0(5) group, and present numerical 
results for the single particle, spin and charge spectra of the SO(5)-symmetric model and identify 
'fingerprints' of SO{5) symmetry in these. In particluar we show that 50(5) symmetry implies a 
'generalized rigid band behavior' of the photoemission spectrum, i.e. spectra for the doped case are 
rigorously identical to spectra for spin-polarized states at half-filling. 



I. INTRODUCTION 



It has recently been proposed Q that the antiferromag- 
netic (AF) and superconducting (SC) phases of the high- 
Tc cuprates are unified by an SO{5) symmetry principle. 
A remarkable degree of continuity between insulating and 
doped phase is indeed supported by a considerable body 
of numerical evidence |^|-^ . Further support for this pro- 
posal came from numerical investigations, which demon- 
strated that the low-energy excitations of physical models 
for the high- Tc materials, i.e. t — J and Hubbard models, 
can be classified in terms of an S'0( 5)— symmetry multi- 
plet structure ^,|^. In these two-dimensional {2D) micro- 
scopic models, the undoped situation - in agreement with 
the experimental situation in the cuprates - corresponds 
to that of a Mott insulator with broken 5*0(3) or spin 
rotational symmetry: long-range AF order is realized. 
The 5*0(5)— symmetry principle then tells us how this 
long-range magnetic order and the accompanying low- 
energy spin excitations are mapped onto the correspond- 
ing off-diagonal long-range SC order and the low-energy 
"Goldstone bosons" (the tt— mode) in the doped situa- 
tion |,|,0. 

However, there exists also a second class of Mott-type 
insulators without long-range AF order, i.e. spin liquids, 
which have a gap to spin excitations. There is growing 
experimental evidence that they are also intimately re- 
lated to the physics of high- Tc compounds: not only do 
these compounds show above the Neel temperature and 
superconducting transition temperature at small dopings 
signs of such a spin gap, but there exist also copper-oxides 
with a Cu02 plane containing line-defects, which result 
in ladder- like arrangements of Ou— atoms P-p^. These 
systems can be described in terms of coupled two-leg lad- 
ders [ p^ , which exhibit a spin gap in the insulating com- 
pound and thus belong to the spin-liquid Mott-insulator 
variety. Also the related "stripe phases" of the 2D Cu02 



planes in the cuprate superconductors have recently re- 
ceived considerable attention jl^Jl^. In these systems, 
the apparent connection between the spin gap and super- 
conductivity must be explained. 

In order to illustrate how the 50(5) theory can, in princi- 
ple, cope with this challenge, an exactly 50(5) invariant 
ladder model has recently been constructed ||l^]. It was 
shown that a 2-leg ladder with entirely local interactions, 
i.e. an on-site interaction \U\ >> t, where t denotes the 
chain hopping, an intra-rung interaction \V\ >> t and 
a magnetic rung-exchange interaction J can have 50(5) 
symmetry if these interactions are related to each other in 
a specific way, i.e. J~4{U+V). The ground-state phase 
diagram of this model was determined in strong coupling 
and, in particular, a regime identified, where the strong- 
coupling ground state is a spin-gap insulator. In addi- 
tion, the spin-gap magnon mode of the Mott insulator 
was shown to evolve continuously into the tt— resonance 
mode of the superconductor. 

However, two key questions remained open, the first be- 
ing the relationship of this 50(5) ladder to the "physical" 
t — J or Hubbard ladders and the second regarding the 
connection to the other variety of Mott insulators, i.e. 
the ones with long-range AF order. 

With regard to the first question, progress was re- 
cently made in the regime of weak-coupling: using the 
weak-coupling renormalization group method, two inde- 
pendent works have recently demonstrated that 
rather generic ladder models at half-filling flows to an 
50(5) symmetric fixed point. 

In the present work, we try to attack both questions 
in the experimentally relevant intermediate and strong- 
coupling regimes. The basic idea is to derive the 50(5) 
concept for a ladder or, more generally, for a spin liq- 
uid with an "RVB vacuum" instead of an "Neel vacuum" 
involving a physically appealing new construction: an 
effective 50(5) invariant low-energy Hamiltonian is con- 



1 



structed in terms of a bond operator representation. In 
strong coupling, this is a well-established concept for spin 
ladders to account quantitatively for the spin ex- 

citations of Heisenberg ladders (up to the physically most 
relevant case of isotropic couplings in the rungs and along 
the ladder). The spin excitations are described in terms 
of triplet fluctuations around the RVB vacuum. The new 
idea is to combine this S'0(3)— symmetric description of 
the spin degrees of freedom , i.e. the triplet excitations, 
with the charge (hole) dynamics. This is accomplished 
by introducing two new Bosons on a rung which stand 
for charge fluctuations, i.e. empty and fourfold occu- 
pied rungs. When combined into a 5-dimensional Boson- 
vector, the dynamics of spin and charge excitations is 
captured by a manifestly SO{5) invariant Hamiltonian. 
It is shown that this formulation in terms of triplet and 
hole fluctuations around an "RVB vacuum" allows for a 
physically transparent demonstration of the corner stone 
in SO{5) theory, i.e. that AF and SC are "two faces of 
one and the same coin". By starting from this "RVB 
vacuum" , which represents the spin liquid state at half- 
filling, we demonstrate that an AF ordered state can be 
generated by forming a coherent state 

|^> ^ e^'^^'^'^irj >, 

which corresponds to z— like triplets condensed into the 
q = TT state. However, in the 5*0(5) theory, the z— like 
triplet with momentum tt and the hole pair with momen- 
tum are components of one and the same 50(5) vector. 
They are rotated into each other by the 50(5) generat- 
ing operator tt. This implies that the above coherent 
state with condensed triplets can be 50(5)— rotated into 
a corresponding coherent state with i^(g = tt) replaced 
by the (hole-) pair creation operator A+(q = 0). This 
state corresponds to hole pairs condensed into the q — 
state, i.e. a superconducting state. In other words: both 
the AF and the SC state can be viewed as a kind of con- 
densation out of the RVB state, or the spin liquid. If 
the so constructed AF state is the actual ground state 
at half-filling, then this physically very appealing 50(5) 
construction yields automatically the ground state in the 
doped situation, i.e. the SC state. 

This construction can also shed light on the second of 
the above questions, namely the interrelation between 
the spin gap Mott-insulator and superconductivity. The 
construction presented here rests on the special geome- 
try of the 2-leg ladder, which suggests a unique "RVB 
vacuum" , from which the AF coherent state can be ob- 
tained. In 2 — Z), such a unique vacuum does, in gen- 
eral, not exist, unless one takes the growing experimen- 
tal evidence for some form of spatial inhomogenities such 
as stripes into account. However, even for a 

translationally invariant state, up to an additional sta- 
tistical average over all possible dimer (singlet) cover- 
ings of the plane, the analogue of the triplet-like exci- 



tation in the ladder can still be generated: here singlet- 
dimers are substituted by triplet-dimers and this excited 
dimer propagates through the 21)— system. Again, the 
50(5) TT- operator converts excited dimers with momen- 
tum (tTjTt) into (i2.2_j^2 —symmetry hole pairs with mo- 
mentum (0,0). A description of spin excitations in 2D 
in terms of Boson-like excited dimers, as for the ladder 
considered here, may thus be a natural starting point for 
clarifying the role of the spin gap (for which the ladder is 
a "toy model" ) and thus the role of the spin- liquid Mott- 
insulators for superconductivity, in general. This will be 
discussed in a forthcoming paper . 
The first question, i.e. the adiabatic connection between 
the 50(5) symmetric model and "physical" ladder mod- 
els, such as the t — J model, is answered numerically in 
the present work. This is achieved by a kind of "Landau 
mapping" of the single-and two-particle excitations of the 
50(5) ladder model onto the corresponding excitations 
of the i — J ladder. 

As mentioned above, the empty-fourfold rung fiuctuation 
involves, in principle, a high energy (of order ~ 2U). In 
the 50(5) description, the Hamiltonian is supplemented 
by a term which contains density and exchange inter- 
actions, and which physically pulls down the empty- 
fourfold fluctuation to be degenerate with the triplet- 
triplet fluctuation. It has been shown in Ref. that 
this can be achieved already for purely local interactions 
V and J within the rungs by choosing 4(J7 + V) = J. 
We demonstrate numerically that, and give physical rea- 
sons why, this at first sight rather unphysical parameter 
condition nevertheless gives in fact an astonishingly good 
mapping of the 50(5) model onto the t— J model for the 
low-energy to versus momentum q dynamics. These find- 
ings strongly support the physical relevance of the 50(5) 
description in physical ladder models, such as Hubbard 
and t — J models, which are of the spin-liquid Mott in- 
sulator variety. 

The paper is organized as follows: Sec. II describes 
the construction of an 50(5) symmetric ladder model 
in terms of bond-triplet and -charge bosons, which are 
unified into a common 50(5) vector (i),.. It also con- 
tains the coherent-state description, where an AF or- 
dered state with staggered magnetization in z— direction 
is generated by a condensation of z— like triplet excita- 
tions. When the triplets are rotated (via the tt- operator) 
onto the corresponding hole pairs (with (1^2 _y2 symme- 
try), a d— SC state results. Sec. HI then contains a gen- 
eral representation of the 50(5) multiplet structure in 
terms of irreducible representations in the Q — Sz plane 
(here Q={Ne — N)/2 with Ne the number of electrons 
and N the sunmber of sites ||^ is the pair density and 
Sz the z-component of total spin). We consider here 
the previously studied multiplets of diamond shape 
(corresponding to even numbers of electrons or holes, i.e. 
electron or hole pairs) as well as multiplets evolving to 
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a square shape. The latter correspond to odd numbers 
of charge carriers and are thus required for building up 
selction rules for photoemission. Sec. IV verifies these 
selelction rules in exact diagonalizations of 2 x 6 5*0(5) 
ladders. Similarly, sec. V numerically illustrates the 
spin and charge dynamics of the exact SO{5) model and 
checks corresponding selection rules for spin- and charge 
correlation functions. Sec. VI discusses the "Landau 
mapping" and adiabatic connection to the physical t — J 
model, and sec. VII gives a short conclusion. 



II. ELEMENTARY EXCITATIONS OF THE 
LADDER 

We begin with the standard Hubbard model on a 2-leg 
ladder, i.e. 
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(1) 



The nearest-neighbor hopping integral along the legs will 
be called t, the one within the rungs t±, and U denotes 
the on-site Coulomb repulsion. We assume the strong 
couphng hmit, U/t,U/t± 1 and start out with the 
case of half-filling and vanishing leg-hopping t. The sys- 
tem then decomposes into an array of Heitler-London 
type dimers and the gound state is simply the product 
of singlets along the rungs of the ladder, with energy 
Eq « —2Nt\/U , where is the number of rungs. 



Singlet 




FIG. 1. Starting from the RVB vacuum (a) the 'virtual' 
hopping process (a) (b) — > (c) leads to the pair creation of 
two triplets (indicated by dashed lines). Further virtual hop- 
ping processes such as (c) (d) ^ (e) enable the propagation 
of the triplets along the ladder 

This state is frequently referred to as the 'rung RVB 
state' and has been the basis of many theoretical works 
on spin ladders |18-2^]. Next, let us consider the effect 
of switching on t. This will produce charge fluctuations, 
and hence exchange processes along the legs. These lead, 
as a first step, to the 'pair creation' of triplets along the 
legs. In subsequent steps, these newly created triplets 



can propagate along the ladder (see Figure ||) . Introduc- 
ing operators which create singlets and triplets : 
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(2) 



we can write the rung-RVB state (which we henceforth 
consider as a kind of 'vacuum'), as 



P) 



N/2 



(3) 



Here, the site indices along a rung, (i,j), have been re- 
placed by a single index n labelling the rung. It has 
been shown by Gopalan et al. that the dynamics of 
the triplets can be mapped exactly onto a system of three 
species of hard-core bond-Bosons on a ID chain. Thereby 
the presence of the Boson t\ on some site means that the 
corresponding rung is in the triplet a state - absence of 
any Boson implies that the rung is in the singlet state. 
The 3 components of the triplet on the n}^ rung form a 
3-dimensional vector 



(4) 



and the Hamiltonian operator governing the triplet-like 
Bosons can be cast into the following manifestly 50(3) 
invariant form 




E 



-t t 

— t t \ . 

''"ri ' '''n+l ''"n+1 ' ''' n ) ■ 



(5) 



Here :: as usually denotes normal ordering and the pa- 
rameters are J—At^/U and J^=At\/U . The terms in (||) 
describe the 'energy of formation' of a triplet (1st term); 
the pair creation of two triplets on neighboring rungs 
(2nd term); the propagation of a triplet (3rd term), the 
simultaneous 'species flop' of two triplets on neighboring 
rungs (4th term); and the exchange of two triplets (5th 
term) 

Let us now consider a new type of fluctuation. Once a 
charge fluctuation has been created one might envisage 
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that the remaining electron in the singly occupied rung 
follows suit, so that the resulting state has one empty 
rung, and a fourfold occupied rung (see Figure ||). Again, 




FIG. 2. A second type of fluctuation which leads to the 
creation of charged rungs: starting from the vaccum, (a) 
(b) (c) leads to the creation of an empty and a fourfold oc- 
cupied rung. Further virtual hopping processes such as (c) — > 
(d) (e) enable the propagation of these charged excitations. 

in subsequent steps the empty and fourfold occupied rung 
can propagate, in much the same way as the triplets did. 
Clearly, the fourfold occupied rung has a very high en- 
ergy K 2U . Following Rcf. |^ however, we now amend 
the Hamiltonian with the goal to 'pull down' this state 
in energy, so as to make the empty-fourfold state degen- 
erate with the triplet-triplet state. This can be achieved 
by adding terms of the form 

^1 = E ^^^^'^^ - " 1) + E -^J' ('^) 

where rii — c| ^Ci^a is the operator of electron den- 
sity on site i, and Si the operator of electron spin. It 
has been shown in Ref. [|l5| that by retaining only a 
density interaction V and exchange constant J within 
the rungs, and choosing U + V=J/A, one can reach a 
situation where the empty-fourfold fluctuation in Figure 
is precisely degenerate with the triplet-triplet fluctu- 
ation in Figure |l|c, and both of them are lower than the 
singly-threefold intermediate state (Figure||3, Figure^D). 
In the following, we call the energy of the singlet Eq, that 
of the triplet /empty/fourfold occupied rung Ei and that 
of a singly/threefold occupied rung £^2- The situation 
we want to reach is Eq < Ei E2- We would like to 
stress that the parameters we need to choose to reach 
this {U + V=J/A) are rather unphysical. Physically, if 
we want J and U to be positive and J not to be too 
large, this requires a large negative rung interaction V. 
However, for the time being we ignore this complication 
and defer a discussion of the case _ffi ^ to the end of 
this section. 

If we restrict the Hilbert space to only rung triplets and 
singlets, the Hamiltonian (ph remains valid also for the 



SO{5) symmetric ladder. The only difference is a change 
of J to J' = / {E2 — El) and a change of the Boson's 
'energy of formation' from ^ Ai = Ei — Eq. Our goal 
is now to simply retain this Hamiltonian, but enlarge the 
3-dimensional vector r into a 5-dimensional vector com- 
prising 2 additional Bosons which represent the empty 
and fourfold occupied rung. We introduce new Bosons h 
and d, whose presence on a rung implies that the rung is 
empty or fourfold occupied: 

\vac) 

(the extra minus sign makes sure that the state corre- 
sponding to dl j can be written as aIjsIj, with A the 
singlet-pairing operator along a rung). Just as two sin- 
glets on neighboring rungs can convert themselves in a 
kind of pair creation process (Figure|l|a Figure|I|c) into 
two triplets, they can also convert themselves into an 
empty rung and a fourfold occupied rung. And, similarly 
as a triplet can exchange itself with a singlet, so can an 
empty or doubly occupied rung. An analogous process 
to the 4th term would be the conversion of two triplets 
on neighboring rungs into an empty and a fourfold oc- 
cupied rung. And finally, an empty rung can exchange 
itself also with a triplet, which corresponds to the 5th 
term, in Eqn.(5). 

The empty rungs and fourfold occupied rungs thus will 
act 'almost' like the triplets - there is, however, a subtle 
and crucial difference. Careful calculation shows that all 
matrix elements for pair creation and propagation of the 
empty and fourfold occupied rungs have the same mag- 
nitude, but opposite sign as for the triplet rungs. The 
ultimate reason is that a spin-flip (which propagates a 
triplet) is accomplished by a back-and-forth motion of 
electrons, and hence picks up a Fermi minus sign, whereas 
the propagation of an empty rung corresponds to a net 
movement of charge in one direction, and hence gets 
no Fermion minus sign. We cope with this by a rung- 
dependent gauge transformation for the charged Bosons 
and introduce: 

tl, = e^--^idi-hi). (8) 

The extra phasefactor e''^" precisely cancels the minus 
sign in the matrix elements of the two new Bosons, so 
that we can now describe the physics of the Bosons by 
the following manifestly 5*0(5) invariant Hamiltonian 

if = Ai^tt 

n 

n n 
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n 

Here t denotes the 5-dimensional vector 
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(10) 



In terms of this 5-diniensional Boson-vector on a single 
rung the 10 root generators of 50(5) take the simple 
form (suppressing the rung index) 



Lab 



tit 



b a/J 



(11) 



which for x < a,b < z reduces to a representation of 
the SO{3) spin operators in the rung basis [^9). Bearing 
in mind that the t are hard-core Bosons, it is straight- 
forward to show that the operators ( |lT| ) obey the 5*0(5) 
angular momentum algebra ||l[: 



[Lab, Led] — ~i{5adLbc — ^acLbd + ^bcLad — SbdLac)- 



Choosing a G x,y, z, the combination 
1 



V2 
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(12) 



(13) 



replaces an a-type triplet with momentum fc = tt by a 
hole pair along the rungs with momentum fc = - this 
operator is therefore nothing but the ladder equivalent 
of the TT-operator in 2D [Q , which replaces a triplet with 
momentum (tt, tt) by a (i^2_j^2 hole pair with momentum 
(0,0). 

We can also express the 50(5) hard-core Boson Hamilto- 
nian (^ directly in terms of a 50(5) quantum nonlinear 
a model Hamiltonian by introducing the 50(5) super- 
spin vector at a given rung as 



and its conjugate momentum as 



(14) 



(15) 



In this representation, the symmetry generator takes the 
more familiar form 



Lab 



XaPb - XbPa 



(16) 



Using these operator identities, the 50(5) hard-core Bo- 
son Hamiltonian (|^) can be expressed as (repeated 50(5) 
indices are summed over!) 



Ai 

4 

4 



(9) i? = ^^i^fc(n) + J'5]a;,(n)a:Jn+l) 

n n 

'Y,Labin)Labin + l) (17) 



This Hamiltonian is quantized using the 50(5) commu- 
tation relations ([121) and 



(18) 



[Lab,Xc] — ~i{SbcXa ~ SacXb), 

together with the hard core constraint 

XaXb — ^ab- 



(19) 



The 50(5) ladder model can be used for a particularly 
simple and transparent demonstration of the key feature 
of the 50(5) theory, namely the one-to-one correspon- 
dence of antiferromagnetism and superconductivity. The 
ground state of the ladder models is actually an RVB 
type of vacuum without AF long-range-order. However, 
for illustrative purposes, let us now construct an AF or- 
dered state (which is in general not a eigenstate of the 
Hamiltonian) by condensing the magnons into the RVB 
ground state. We can obviously express the operator of 
staggered magnetization in z-direction as 



AL 



e*™(/^„(Ti)-P„(iT)), 



where e.g. -P„(ti) projects onto states where the n*-^ rung 
has the configuration t|. It is now easy to see that 



whence 



{Pn{n)-Pn{i]) )4 = 4., 

{PnViD-PniiX) )tla^Sa,. 4: 

= tl{q = n)+t^iq = 7r)]. 



(20) 



If we now form the coherent state 



^\^/Ntliq=■^) 



which corresponds to z-like triplets condensed into the 
fc=7r state, and treat the t as ordinary Bosons, we obtain 

(^'aIM.I^'a) = V2XN. 

If the hard-core constraint is taken into account rigor- 
ously, the only change is an extra correction factor of 
1/(1 -I- A^) on the r.h.s., see the Appendix. 
This calculation shows that by starting from an 'RVB 
vacuum', an antiferromagnetically ordered state with 
A/5 in z-direction can be generated by condensing z-like 
triplet-excitations into the k—n state. At this point we 
can invoke the 50(5) symmetry of the model, which tells 
us that since the z-like triplet with momentum tt and the 



5 



hole pair with momentum are two different components 
of a 5-vector (the difference in momenta is precisely ab- 
sorbed by the gauge transformation (^) we needed to 
make the signs of hopping integrals consistent), they are 
dynamically indistinguishable. This means that the AF 
state, with condensed triplets, can be S'0(5)-rotated into 
a state with condensed hole pairs. It follows that if the 
antiferromagnetic state were the ground state at half- 
filling (which is the case for 2D materials) we can replace 
all z-like triplets by hole pairs with momentum and by 
50(5) symmetry automatically obtain the ground state 
in the doped case. The latter then consists of hole-pairs 
along the rungs condensed into the k = states and thus 
is necessarily superconducting. In other words: both the 
antiferromagnetic and the superconducting state may be 
viewed as some kind of condensate 'on top of the rung- 
RVB state. 5*0(5) symmetry then simply implies that 
the condensed objects are combined into a single vector, 
whence the unification of antifcrromagnetism and super- 
conductivity follows in a most natural way. 
The above derivation makes sense only in a strong cou- 
pling limit, where a discussion starting out from rung- 
singlets makes sense. One might expect, however, that 
similar considerations will apply also for cases with a 
weak coupling within the rungs Jl6| , pTt . 
To conclude this section, we discuss what will happen 
if we switch off the correction terms (^, which were in- 
troduced so as to enforce exact 50(5) symmetry. The 
crucial question is whether above considerations will re- 
tain some validity or break down alltogether. Let us 
consider the 1 — 5 plane of the 5-dimensional space, i.e. 
the 'charge-like subspace'. In the fully 50(5)-symmetric 
model the vector t can be rotated freely in this plane. 
The total charge is related to the angular momentum in 
the 1 — 5 plane, since Q = L15 = Xip^ — x^pi. For a sys- 
tem with a large charge gap compared with the spin gap, 
one can apply a chemical potential to lower the energy 
of the empty rungs hl^\0) at the cost of further increas- 
ing the energy of the fourfold occupied rungs dl^\0). The 
chemical potential acts like a magnetic field normal to 
the 1 — 5 plane, and selects a particular sense of the rota- 
tion in the 1 — 5 plane. Therefore, the "free rotor" in the 
exact 50(5) symmetric models becomes a "chiral rotor" 
in the more realistic Hubbard or t — J models. 
In terms of the Boson model, switching off Hi will lead 
to different energies for the different Boson species, i.e. a 
replacement of the form 

n n 

where eo corresponds to the binding energy of a hole pair 
(relative to the rung singlet) and £2 ~ 2J7. Here we are 
neglecting the modification of the off-diagonal matrix el- 
ements cx J'. While 62 is a huge energy, let us assume 
that we have a state with Q < and treat this term in 
perturbation theory. The crucial point is now, that for 



total electron density < 1 the fourfold occupied rung dj^ 
is admixed only as a quantum fluctuation, so that we 
expect for the change in energy 

In the limit J' <C i?i — i?o the main effect of switching 
to the physical Hubbard ov t — J model in the hole-doped 
subspace thus is adding a chemical potential-like term 
(eo ~ Ai)Q to the Hamiltonian. Apart from that, hav- 
ing 62 2U for states with Q < merely corresponds 
to projecting out quantum fluctuations involving charged 
Bosons. This may be expected to have only a minor ef- 
fect. Of course states with Q > 1 will have their energies 
shifted by Qe2, so that 50(5) rotations into this sector 
of the Hilbert space become forbidden. As long as we 
restrict ourselves to hole-doped states, however, one may 
hope that 50(5) symmetry remains approximately valid. 
Below we will present an explicit numerical check of this 
more qualitative discussion. 

III. CLASSIFICATION OF STATES ACCORDING 
TO SO(5) MULTIPLETS 

In this section we shall briefly review the 50(5) group 
theory and classify all states in the ladder model accord- 
ing to irreducible representations (irreps) of the 50(5) 
Lie algebra. We also discuss various 50(5) selection rules 
which shall be used in the next two sections. 
In reference Q , three of us showed how low lying bosonic 
states in the t — J model can be classified by the fully 
symmetric tensor multiplets of the 50(5) Lie algebra. 
The full multiplet structure of the 50(5) group is much 
richer. The general irreps of an 50(5)-symmetric model 
are described by two integers 

{p,q),P>q>0 (21) 

with dimension 

D = il + q){l+p-q) (l + ^p) (1 + J(P + '?)) , 

(22) 

and Casimir ||2|] 

C = ^p^ + ^q^ + q + 2p. (23) 

50(5) is a rank two algebra so we choose the charge Q 
and the z-component of the Spin 5^ as members of the 
Cartan subalgebra of mutually commuting generators. 
This allows us to draw the different irreps in the Q — Sz 
plane. The (p,p) irreps have the familiar ||^ diamond 
shape and corresponds to the fully symmetric traceless 
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tensor irreps identified in . The {p, 0) irreps form a 
square in the Q ~ Sz plane and are the spinor states 
found in systems with odd number of electrons. For ex- 
ample, the quintet Ei manifold with 3 magnons and 2 
pair states carry the label (1,1), while the quartet E2 
manifold with one or three electrons per rung carries the 
label (1,0). The intermediate multiplets {p,p > q > 0) 
are the evolution from the square to the diamond (Figure 
All states of the exact 5*0(5) symmetric models 
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FIG. 3. The first few irreps of SO{5) displayed on the 
Q — Sz-p\a,ne. The {p,p) multiplets have the familiar diamond 
shape; with {p, q < p) the multiplets evolve to a square shape. 
The number of dots in each multiplet is not the full dimension 
of the 50(5)-irrep since there are additional degeneracies in 
the Sz-direction. 

constructed in |l^ can be classified into this 5*0(5) multi- 
plet structure. We have verified this numerically by com- 
puting the energies and expectation values of the Casimir 
operator for all eigenstates of a 2 x 4 system. By scan- 
ning all sectors of different hole number and z-spin and 
collecting states of equal (within computer accuracy, i.e. 
10^^^) energy and Casimir charge we have verified that 
from the ground states up to the highest excited states all 
eigenstates of the system can be classified into the mul- 
tiplets shown in Figure || (and more complicated ones). 
For the 2x6 system, where a full digonalization is not 
feasible anymore, we have verified this for the low energy 
states obtainable by Lanczos diagonalization. 
Besides classifying eigenstates into 5*0(5) multiplet 
structures, 50(5) symmetry also gives powerful selec- 
tion rules on the possible tranition processes and relates 
various matrix elements through the Wigner-Eckart the- 
orem. In the following two sections, we are interested in 
the photoemission, charge and spin spectra of the 50(5) 
symmetric ladder model. The perturbing operator in the 
photoemission process is a single electron operator, which 
transforms according to the 4 dimensional (1,0) irreps 
under 5*0(5). In this work, we shall consider the initial 
and final states of the form (p,p), corresponding to ex- 
cited magnon and pair states. The possible final states 
generated by the perturbing operator can be obtained by 
decomposing the product representation into irreducible 



components, and it given by 

{l,0)®{p,p) = {p+l,p) + ip,p-l) (24) 

Similarly, the perturbing operator in the neutron scatter- 
ing or Josephson tunneling process transform according 
to the 5 dimensional (1, 1) irreps under 50(5). The pos- 
sible final states are given by the following decomposition 
rules: 

(l,l)(g)(p,p) = (p+l,p+l) 

+ ip+l,p-l) + {p-l,p-l) (25) 

The photoemission selection rules (except for 5 — > 
5 ± i) can be easily visualized by superimposing the 
involved 5*0(5)-multiplets (Figures ^ and ||). Figure ^ 
shows that there are four possible photoemission and in- 
verse photoemission transitions for the Sz = state at 
half- filling which belongs to the (1,1) multiplet. We 
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FIG. 4. Visualization of the selection rules for photoemis- 
sion and inverse photoemission for 5*0(5) symmetric states on 
the exampfe of a transition originating from a haff fiffed state 
with 5*z = in the (1, 1) multipfet. SO{5) selection rules only 
alfow transitions to the (1,0) and the (2, 1) multiplets (black 
discs), which are superimposed on the initial (1, 1) multiplet 
(circles). In this example, the photoemission and inverse pho- 
toemisson process remove/inject a spindown-electron so the 
photoemission transition is visualized by an arrow pointing 
south-east and the inverse photoemission arrow points to- 
wards north-west. 

therefore expect four distinct bands in the spectrum. If, 
on the other hand, one takes a spin-polarized or doped 
state of the same multiplet as the initial state for pho- 
toemssion and inverse photoemission we expect to see 
only three bands since there are only three possible tran- 
sitions according to the selection rules (Figure Fur- 
thermore, we expect the very same spectrum for both 
the polarized and the doped case since the initial states 
belong to the same initial multiplet as well as all tar- 
get states belong to the same target multiplets and are 
therefore degenerate. 
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FIG. 5. Photoemission and inverse photoemission originat- 
ing from a spin-polarized half-filled state {Q = 0,Sz = +1) 
and from a doped state (Q = —1,5*2 ~ 0). Both states are 
members of the (1, 1) multiplet like the initial state in Figure 
4. This Figure shows that there is no allowed photoemission 
transition to the (1,0) multiplet for either of the two states. 



IV. SINGLE PARTICLE SPECTRA 

We now want to see the practical application of the 
selection rules discussed in the preceeding section. To 
that end we study the photoemission (PES) and inverse 
photoemission (IPES) spectrum, which we define as 



ApEsik,^,) = i3(0|4 — — i 

TT ^ uj + H — en — lO^ 



AiPEsik,u;) = -3(0|cfc 



1 



l,|o>- 



Here |0) is a suitably chosen initial state and ep its en- 
ergy. In the following, we will label different states as 
™i^Q, where m denotes the standard spin multiplicity, K 
the total momentum and Q the charge quantum number. 
For finite clusters of our 50(5) symmetric ladder model 
this is readily obtained numerically by Lanczos diago- 
nalization ||2^. The 50(5) multiplets discussed above 
are easily identified, because the energies of states be- 
longing to one multiplet are degenerate (i.e. identical to 
computer accuracy). This allows to study the evolution 
of the single-particle spectral function as we pass from 
one multiplet (p, q) to the other, and within one multi- 
plet through the different doping levels. To begin with. 
Figure || shows the single-particle spectrum for the half- 
filled ground state ^(0, 0)o, which actually forms a one- 
dimensional (0, 0) multiplet. Final states can only belong 
to the 4-dimensional (1,0) representation, see Figure ||. 
Despite the fact that we are using very strong interaction 
parameters, there is just one single electron-like band in 
PES, whose dispersion closely follows the noninteracting 
dispersion. The center of gravity of this band is given by 
the energy difference between a rung-singlet and a singly 
occupied rung, i.e. Eq - E2 = -7U/2 - 3V, i.e. -10 
with our present parameter values. We proceed to the 
spectrum of the ■^(tt, 7r)o state, with 5z=0 - this state 




FIG. 6. Comparison of the single particle spectral function 
for the half-filled ground state (dotted line) and the half-filled 
^(7r,7r)o state with Sz=0 (full line). The Fermi energy, de- 
fined as the average of first ionization and affinity energy, is 
taken as the zero of energy. The parameter values are t = 1, 
t± = l,U = 8,V = -6, and J = 8. 



belongs to the fivefold degenerate {p = l,q = 1) multi- 
plet. From Figure ^ we expect final states belonging to 
both the (2, 1) and the (1, 0) representation. It appears as 
if the bands seen in the ground state spectra remain prac- 
tically unchanged - as a new feature, howewer, there ap- 
pear some weak 'sidebands' close to fi. They can be seen 
both in photoemission and in inverse photoemisson, so 
that we obviously have precisely the 4 bands expected on 
the basis of our discussion of selection rules given above 
(see Figure ^). The physical interpretation is straight- 
forward. The ^(7r,7r)o state has a single triplet-like Bo- 
son. The 'main bands', which are similar to those in the 
ground state spectra, correspond to final states where the 
photohole is generated in a singlet-rung. They therefore 
contain the original triplet, plus a singly occupied rung 
which propagates through the ladder and carries the en- 
tire momentum transfer. These states therefore belong 
to the (2, 1) representation. The center of gravity of this 
band is again at — E2. The two excitations (triplet 
and hole) now can scatter from each other, whence the 
respective band becomes broadened. There is however 
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also a second process, namely the photohole can be cre- 
ated in the rung occupied by the triplet. This creates fi- 
nal states containing only a singly occupied rung, which 
must therefore belong to the (1,0) representation. As 
the intensity of this second process is proportional to the 
triplet density, our interpretation of the sidebands can 
be checked numerically by comparing the mean weight 
of the sidebands for one-triplet-states on ladders with 
different numbers of rungs. On 4-, 6-, and 8-rung lad- 
ders (which corresponds to triplet densities of 1/4, 1/6, 
and 1/8) and with our parameter values, we found mean 
weights of 0.276, 0.183, and 0.137 - a convincing proof of 
the interpretation given above. The center of gravity 
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FIG. 7. Comparison of the electron removal spectrum 
for the two- hole ground state "^(0, 0)_i (dotted line) and the 
half-filled ^(7r,7r)o state with & = 1 (full line). The removed 
electron has |-spin. The spectra for the half-filled state have 
been ofi'set in j/-direction so as to faciliate the comparison. 
Parameter values are as in Figure 6. 

of the resulting band is at the energy difference between 
a rung triplet and a singly occupied rung, Ei — E2 = 
U /2 + V, i.e. —2 with our parameters. In this type of 
process the photohole has to absorb also the momentum 
of the triplet, (tt, tt), so that the dispersion of this band 
is precisely the same as that of the 'main band' seen 
for the ^(0,0)0 ground state, but shifted by (7r,7r) (note 
however the rigid shift of the band by Ei — Eq due to 
the difference in initial state energy). Literally the same 
arguments hold for the inverse photoemisison part, and 



inspection of Figure |6| shows that all of these features are 
observed. 

So far we have seen that spin excitations in the ground 
state generate additional sidebands in the single parti- 
cle spectrum. At this point, however, we can invoke the 
exact SO{5) symmetry of the model, which tells us that 
spin polarization and hole/electron doping are equiva- 
lent, in that the empty rung is the 'SO{5) partner' of the 
triplet rung. Consequently, Figure ^ compares the sin- 
gle particle spectrum of the ^(0, 0)_i and '^(tt, 7r)o states. 
Both states belong to the same (1,1) representation - un- 
ambiguous evidence is provided by the fact that their en- 
ergies agree to computer accuracy. In Figure we have 
chosen the S'x=l member of the triplet, and study the 
annihilation of a |-electron. ^From Figure ^ we expect 
that this simulation probes exclusively the (2, 1) repre- 
sentation. Since the final states for these two different 
processes belong to the same 50(5) multiplet, we expect 
their amplitudes to be directly related. More precisely 
we have \^{tt, 7t)o) = vr^l^O, 0)_i), with tt^ = 7r| -I- in^, 
whence: 

e{TT,n)o,S, = l|4i— 4 Cki\H^,^)o,S. = 1) 

* uj + u — eo 

= (i(0,0)„i|^-,4^— CH^l\HO,0)-i) 

= {\0,0)_,\ci^—-l [^+,4] cfcji(0,0)_i), 

= C(0,0)„i|4— c;,Ji(0,0)_i), 

Where we have used the fact that the 7r+ operator com- 
mutes with Cfc| and the Hamiltonian. Furthermore, it 
annihilates the |^(0,0)_i) state. Since they belong to 
the same S0{5) multiplet, the energies of the |^(7r, 7r)o) 
and the |°(0,0)_i) state are identical, i.e. cq. There- 
fore, as we indeed see in Figure 0, these two correlation 
functions are completely identical (within computer ac- 
curacy in the simulation). The physical reason is that 
the photohole cannot be created in the triplet rung (be- 
cause the latter contains two "f-electrons). In other words, 
when 'seen through the eyes of the photohole-operator' 
the triplet rung looks like it were 'empty', i.e. occupied 
by holes. 50(5) symmetry thus implies that the single 
particle spectra in the doped ground states actually are 
identical to those of certain states at half-filling. SO (5) 
thus presents a very obvious rationalization for the rigid- 
band behavior observed numerically in the single particle 
spectra of the 2D t — J model |^,|| . 
Next, we return to the photoemission spectrum for the 
5^=0 member of the ■^(tt, 7r)o state. In this case, the 
photohole can be generated in the triplet rung, and, as 
discussed above, the sidebands near fj, appear. The corre- 
sponding final states have only one singly occupied rung 
remaining in the system. Then, precisely the same kind 
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Remarkably enough, the photoemission spectrum Apes 
of the half-filled triplet Sz = state is related to the 
photoemission spectrum A'p^g and inverse photoemis- 
sion spectrum A'jp^g of the two hole ground state. In- 
spection of Figure ^ shows the three spectra in question 
demonstrates that the photoemission spectrum for the 
doped ground state (full line) indeed is simply the mean 
of the two spectra computed for the doped case. 
All in all, the preceeding discussion has shown that 
S0{5) symmetry leads to at first sight unexpected be- 
havior in the single-particle spectral function: since spin- 
polarization and hole doping are equivalent under SO{5), 
the PES spectra of half-filled but spin-polarized states are 
identical to those of doped states. This may be the key 
to understand the rigid-band behavior observed (2|j|] in 
the 2D t — J model. Moreover, sidebands which appear 
in PES on spin-excited states at half-filling and which 
originate from processes where a spin excitation is anni- 
hilated, reappear in IPES on hole-doped states. Again, 
this was found previously also for the 2D t—J model [24[. 
We note that the spectroscopies shown in Figures [£ TQ; 
and H have been performed for the actual 2D t—J model 
as well and have produced results in strong support of 
50(5) M. 
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FIG. 8. Comparison of the electron removal spectrum for 
the two-hole ground state A'p^g{k,ij) (dotted line), the 'in- 
verted' electron addition spectrum for the two-hole ground 
state A'jp^g{—k + Q, —u>) (dashed line) and the electron re- 
moval spectrum ApEs{k,uj) for the ^(7r,7r)o state with 5*2=0 
(full line). Parameter values are as in Figure 6. 

of state results if a photoelectron is created in an empty 
rung. The only difference is, that in the former case the 
photoholc has to absorb the momentum of the triplet i.e. 
(tt, tt). We can conclude that the sideband seen in PES at 
half-filling must disappear in PES for the two-hole case, 
but reappear in IPES, shifted by (tt, tt). This can also be 
seen explicitly by the following identity: 

Cin, 7r)o, 5, = 0\cl^—l—Ckif{7T, 7r)o, & = 0) 



H 



C{0,Q)-i\nzcl 



1 



= JC(0,0)-i|4r 



H 
1 



£0 



-Cfcx7rt|i(0,0)_ 



H 



-Hi(i(0,0)_i|c_fe+Qt; 



- eo 
1 



-cLfc+QTl'(0'0)-i>' 



H-eo 

where the second term in the last equation arises from 
the nonvanishing commutator between tTz and Cki- This 
identity implies that 

ApEs{k,uj) = i [A'pEsik^uj) 



(26) 



V. SPIN AND CHARGE DYNAMICS 



Our rung-Boson model gives a description of the 
physics in terms of spinful and charged 2-electron ex- 
citations, and thus should be directly applicable to a dis- 
cussion of the spin and charge dynamics. We introduce 
the spin correlation function (SCF) and the density cor- 
relation function (DCF) as 



D,{k,Lo) = -%Q\S\ ^Sm, 

TT W — il + £0 — 

D,{k,oj) = -5(0|n_fc 7^"fc|0), 



where and n denote the operator of z-spin and elec- 
tron density. As a first step we need representations of 
these operators in terms of the rung-Boson operators t. 
Following Ref. [|9|, the spin operators for the v}^ rung 
and the two possible momenta in y-direction are: 



5(71,0) =XT X T„ 

S{n,TT) = — (rt +-r„). 



(27) 
(28) 



The spin operator with transverse momentum thus ac- 
tually corresponds to a two-particle spectrum, and thus 
will take (in an infinite system) the form of a continuum, 
whereas the spin operator with transverse momentum tt 
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is a single particle-like spectrum. This difference can in- 
deed be seen in the numerical spin correlation function 
(Ref. (19|). 

The 'S0{5) partner' of the operator S{n,TT) therefore is 



1 



(29) 



which is related to the creation and annihilation of a 
singlet electron pair. On the other hand, the operator 
of electron density with transverse momentum is the 
'50(5) partner' of 5(fc,0): 



(30) 



and t'^ 5 were defined in (|g) . Comparing ( 30 ) 



where t^^ a i 

( p7| ) to (|l]) we note that the density and spin opera- 
tors for transverse momentum transfer are actually root 
generators of SO{5) - it follows that we can invoke the 
Wigner-Eckart-theorem (or intuitive arguments) to show 
that the respective spin and density spectra for certain 
pairs of states are identical. Below we will present vari- 
ous examples. 

As for the density operator with transverse momentum 
transfer tt, this is much more complicated because this 
operator annihilates any state corresponding to either 
charged or spinful Bosons. This correlation function 
therefore will be determined by the dynamics of the 
singly and threefold occupied rung, which in our strong 
coupling limit is admixed only virtually. We therefore 
expect that this correlation function has only weak fea- 
tures at relatively high energies. 

After these preliminaries we proceed to a discussion of 
the numerical spectra. With regards to the SCF it is im- 
portant to keep in mind that all following spectra have 
been computed using the z-component of the spin oper- 
ator. 

The half-filled ground state does not contain any Bosonic 
excitation at all. We can obtain a nonvanishing SCF only 
by coupling to the weakly admixed singly and threefold 
occupied rung, and therefore expect an extremely weak 
signal for kT=Q- Final states for the SCF with kx—TT 
are in the (1,0) representation, and the dominant peak 
in the respective spectra therefore gives the dispersion of 
the triplet Boson (see Figure P). Choosing now the low- 
est of these (1, 1) states, the ^(tt, 7r)o, as the initial state 
for the calculation of the SCF with kT=0, the triplet Bo- 
son with momentum fc^vr could in principle be scattered 
into a state with momentum fc + tt, where k is the mo- 
mentum transfer along the ladder. Since we are using 
the Sz operator, however, this is not possible (see the 
representation (p7|)) if we use the Sz—0 member of the 
'^(7r,7r)o state. We therefore expect a very low intensity 
spectrum which again is due to virtually admixed singly 
and three-fold occupied rungs. Moreover, since the z-like 
triplet acts precisely like a hole pair, and since the hole 
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FIG. 9. Comparison of the spin correlation spectrum for 
the half-filled ground state (dotted line) and the "^(vr, 7r)o state 
with 5*2=0 (full line). Parameter values are as in Figure 6. 



MO 




FIG. 10. Comparison of the spin correlation spectrum for 
the two-hole ground state (dotted line) and the half-filled 
^(7r,7r)o state with Sz = l (full line). The spectra for the 
two-hole state have been offset in y-direction to facilitate the 
comparison. Parameter values are as in Figure 6. 
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pair cannot be excited by the spin operator either, the 
kT=0 spectra of the •^(7r,7r)o state with Sz=Q and the 
^(0,0)_i ground state must be identical by 5*0(5) sym- 
metry (which indeed they are, compare Figure p] ). The 
situation changes completely if we use the Sz — l member 
of the •^(7r,7r)o triplet as initial state (see Figure p^): this 
state contains a mixture of the x-like and y-like triplet, 
which both can be excited by the .Sz-operator with kT=0- 
We will therefore see a strong peak, with essentially the 
same dispersion as the dominant peak in the fcT=7r spec- 
tra in Figure ^, but shifted by Ak=iT in momentum and 
by w _Ei — Eq in energy. On the other hand acting with 
the z-like spin operator with kT='!r creates or annihilates 
a z-like Boson. Then, in both states, "^(tt, 7r)o with 5*2=1 
and ^(0, 0)_i, annihilation of a 2;-like Boson is not pos- 
sible. Creating a z-like Boson, we will generate 50(5) 
equivalent states, because the created z-like Boson inter- 
acts with the one already present in the system (a;-like, 
y-like or hole- like) in identical ways - the kT=T^ spectra 
for these states thus must be identical, which indeed they 
are (compare Figure |l^). Comparing Figure |l0| with Fig- 
ure we see that the magnon mode with (tt, tt) at half- 
filling evolves continuously into the tt resonance mode of 
the two hole state, in accordence with the analysis of ref- 
erence [ p^ . Away from (7r,7r), a lower energy structure 
emerges in the SCF of the two hole state. 
On the other hand, when we switch to '^(7r,7r)o with 
Sz—0, the parts corresponding to a creation of the z-like 
Boson are slightly different from the spectra in Figure |l^ 
because the final state now contains two Bosons of equal 
(z-like) species, which interact differently as compared to 
unequal species. 

We proceed to a discussion of the DCF. In the remainder 
of this section we will refer only to spectra with kx—O. It 
follows by 50(5) symmetry and (|3), that the den- 
sity operator with kT=0 takes a completely analogous 
form as the spin operator. When using an 5*0(5) singlet 
as initial state, both operators thus must give the same 
spectrum. This can be seen by comparing the SCF and 
DCF for the half-filled ground state, compare Figure ^ 
and |l2|. They are also identical when acting onto states 
which contain only z-like Bosons - which can be excited 
neither by the 52 or the density operator. This explains 
the identity of the SCF and DCF for the "^(tt, tt) state 
with 52=0, compare Figure |9| and Figure 
On the other hand, since the hole pair behaves like a 
triplet, the DCF for the two- hole ground state must be 
identical to the SCF for the ^(tt, 7r)o state with 5*2 = 1, 
compare Figures |l^ and |l2|, right panel. 
Summarizing this section, we have seen that 50(5) sym- 
metry enforces identity relationships between spin and 
charge correlation functions for different initial states, 
which differ by their hole number. Again this implies 
strong continuity with hole doping. While such continu- 
ity is certain to be present in any system with (approxi- 
mate) 5*0(5) symmetry, it should be noted that the 
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FIG. 11. Comparison of the spin correlation spectrum for 
the two-hole ground state (dotted line) and the ^(7r,7r)o state 
with 5*z=0 (full line). The spectra for the two-hole state have 
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FIG. 12. Comparison of the density correlation spectrum 
for the half-filled ground state (dotted line, left), the ^(tt, 7r)o 
state with 5*z=0 (full line, left), and the two-hole ground state 
(right). Parameter values are as in Figure 6. 
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identity relationships in the present case depend very 
much on the representations (|2^) and (^) of the spin 
and density operators in terms of the rung Bosons. While 
these are easily transferable to ladder-like systems, anal- 
ogous identities for a fully planar system may be quite 
different and more restricted. 



VI. ADIABATIC CONNECTION TO THE T-J 
MODEL 

So far, we have studied the spectra for the ideal SO{5) 
symmetric model. In this section we want to investigate 
the effect of discarding the correction terms Hi which 
enforced exact 5*0(5) symmetry. More precisely our goal 
is to find, in the spirit of Landau's Fermi liquid theory, a 
'mapping' between the excitation spectrum of the exactly 
50(5) symmetric but in principle unphysical model, and 
that of a 'physical' t — J ladder. The latter, being the 
[7— >cxD limit of the Hubbard model, incorporates the con- 
straint of no double occupancy which, based on erro- 
neous mean-field calculations, has recently been argued 
to break 50(5) symmetry |2^]. We compare the single 
particle spectra as well as the spin correlation (i.e. a two 
particle spectrum) of the 50(5) symmetric ladder and 
the t — J model for a 6-rung ladder. Thereby we adjust 
the parameters of the 50(5) symmetric ladder so as to 
obtain an optimal fit to the t — J model, for which we fix 
t±/t—Q.b and J/t—Q.5, J±/t=2.0. We are using a rela- 
tively large value of the exchange along the rungs, J± - 
this is by no means adequate to describe actual materials, 
but our main goal here is to investigate the effect of the 
constraint of no double occupancies in the t — J model. 
Note that the t — J model has a 'Hubbard gap' which is 
effectively infinite, whereas in the 50(5) symmetric lad- 
der hole-doped to electron-doped states are degenerate. 
The models thus appear at first sight entirely unrelated, 
but our goal is to check our above conjecture that as long 
as we restrict ourselves to the hole doped sector, the two 
models still can well have essentially the same low-energy 
dynamics. 

Figure [l^ and Figure |lj show the photoemission spec- 
tra originating from the half-filled ground state ^(0, 0)o, 
and ^(tt, 7r)o with 5^=0, The overall agreement of the 
two model's spectra was obtained choosing U/t=A, 
V/t=-3.25, J/t=3, and tj_/t^2/5 for the exactly 
50(5) symmetric model, followed by an energy rescal- 
ing by a factor of 0.625 (which amounts to setting 
^so(5)— 0-625tt_j). The resulting mapping is excellent: 
the t ~ J model's number of peaks, their weight distri- 
bution and dispersion are accurately obtained also in the 
exactly 50(5) symmetric model. In particular, the ap- 
pearance of the 'sidebands' in the spectra for spin polar- 
ized ground states can be nicely seen for both models. 
Regarding the spin correlation for the same set of pa- 
rameters (Figure 15), there is still good agreement as to 



the 



k = 




FIG. 13. Comparison of the photoemission spectra from 
the half-filled ground state of the t — J model (dotted line) 
and the exact 50(5) model (full line). The t — J parame- 
ters are t±/t=G.b, J/t—0.5, J±/t—2 , the SO{5) parameters 
1 1 /t=2/3. U/t=A. Wt=3.25. J/t=3. with tc:n,K-.=0.625tt-j. 



MO 






FIG. 14. Comparison of the photoemission spectra from 
the half-filled ^{n, n) state with Sz=0 of the t — J model (dot- 
ted line) and the exact 50(5) model (fuU line). Parameters 
are as in Figure 13. 
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number and the weight of the peaks and the character 
of the dispersion. The absolute size of the dispersion is 
somewhat smaller in the exact S0{5) model but the value 
of the 'spin gap' at (tt, tt) is reproduced quite accurately. 
Alltogether, the mapping between the 'unphysical' but 
exactly SO(5) symmetric model and the physically bet- 
ter founded t — J model is reasonably good, which clearly 
supports the physical relevance of SO{5) symmetry. We 
also note that the calculations are 'as close as possible' 
to half-filling, where the impact of the Hilbert-space pro- 
jection still is presumably the largest in the t — J model - 
yet the agreement of the low-energy physics is obviously 
quite good. One therefore may hope that similar map- 
pings can be carried out also for more realistic parameter 
values oi t — J ladders, so that for f — J and Hubbard 
ladders the SO{5) symmetric model may in fact be the 
generic effective Hamiltonian. This is further supported 
by numerical studies in the case of large U^—V, J=0 
|p7| recently carried out by Duffy, Haas and Kim, and 
renormalizatiion group calculations |p^ , |l7| for the weak 
coupling case. 



. t-J-Model 
_ exact S0(5) 




40 
(o/t 



FIG. 15. Comparison of the spin correlation spectra for the 
half-filled ground state of the t — J model (dotted line) and 
the exact SO{5) model (full line). Parameters are as in Figure 



VII. CONCLUSION 

In summary, we have studied the dynamics of an ex- 
actly 5*0(5) symmetric ladder model. We have shown 
that in the strong-coupling limit the model reduces to 



the simplest possible SO (5) symmetric generalization of 
a Bosonic Model derived previously by Gopalan et al. 
I p^ for the spin excitations of Heisenberg ladders. In this 
limit the 'vacuum' of the theory is the rung-RVB state, 
and the elementary excitations of the 5*0(5) symmetric 
ladder correspond to uncharged triplet-like and charged 
singlet-like Bosons. 50(5) symmetry then simply states 
that the spin-like and charge-like Bosons are dynamically 
equivalent, in the same sense as e.g. proton and neutron 
are considered dynamically equivalent in the isospin the- 
ory of nuclear physics. In the strong coupling limit, the 
low energy physics can be mapped onto a model of five 
species of hard-core bosons. The effective non-linear a 
model description j^] can be systematically derived in 
this limit. 

Because the ground state of the ladder system is a Mott 
insulator with 'RVB type" of singlet vacuum, it can 
serve as a reference state to consider the condensation 
of magnons and charged bosons on the equal footing, 
thereby revealing the precise analogy between antiferro- 
magnetism and superconductivity. In a 2D system, no 
such reference state is found to be the ground state for 
any reasonable Hamiltonian in the infinite system. How- 
ever, the ground state of any finite cluster is a total spin 
singlet, and recent numerical calculations demon- 
strate that the ground state of the t — J oi Hubbard 
model is also approximately a 5*0(5) singlet. The low 
energy excitations of the 2D cluster are magnons and 
hole pairs, and share many similarities to the properties 
of the 50(5) ladder system found in this work. 
As shown in the present work, 50(5) symmetry has pro- 
found implications for the dynamical correlation func- 
tions, most notably the single particle spectrum. Specific 
predictions of 50(5), like a 'generalized rigid band be- 
haviour' I^J^ and the appearance of sidebands in the in- 
verse photoemission spectrum may indeed have been 
observed long ago in the actual 2D t — J model. Mo- 
tivated by the present theory we have carried out more 
detailed spectroscopies on the 2D model and obtained 
results in strong support of 50(5) 
Finally we (and other authors 
strated that despite the at first sight rather unphysical 
parameter values of the 50(5) symmetric model, a 'Lan- 
dau mapping' to the more realistic t — J model is feasible. 
This may suggest that the 5*0(5) symmetric ladder is in 
fact the generic effective Hamiltonian for two-leg ladder 
systems. 

Acknowledgements: Two authors (W.H. and S.C.Z.) have 
benefitted from many enlightening discussions with D.J. 
Scalapino. SCZ is supported in part by the NSF under 
grant numbers DMR94-00372 and DMR95-22915. A.D., 
M.G.Z. and W.H. acknowledge support from DFN con- 
tract TK598-VA/D3 and are grateful to the University of 
Stuttgart and Miinchen Supercomputing Centers. R. E. 
most gratefully acknowledges the kind hospitality at the 
National Center for Theoretical Studies, National Tsing 




) have demon- 



14 



Hua University, Taiwan. 

VIII. APPENDIX 

We consider the coherent state 

J_„AvWt,(7r)t. 



vac) 



where n is the normahzation factor, and t2(7r)^ the cre- 
ation operator for a true hard-core Boson. Let us assume 
that the exponential has been expanded, and ask for the 
total norm, rt^, of all states of v*^ order in A (i.e. all 
states containing v Bosons). Each of these terms has a 
prefactor of y /v\ from the expansion of the exponential. 
There is also a factor of ±1, depending on how many 
Bosons are on rungs with odd numbers - this will disap- 
pear upon squaring the prefactor and we therefore disre- 
gard it. Moreover, the factor of l/v\ is cancelled because 
each Boson configuration is generated in v\ different ways 
when expanding (tz{'K)'^Y ■ The number of different Bo- 
son configurations (which are mutually orthogonal!) with 
V Bosons is Nl/{N — lyy.i'l, whence 



A 



2iy 



and summing over u we obtain 

n= (l + A^)^. 

Next, we have 
whence 

{^x\Ms\^x) = -^9Alog(n) 
V2\N 



1-f A2' 



q.e.d. 
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